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Summary 

The dynamic response of a fiber optic Bragg grating to mechanical vibrations is examined both 
theoretically and experimentally. The theoretical expressions describing the consequences of changes in 
the grating’s reflection spectrum are derived for partially coherent beams in an interferometer. The 
analysis is given in terms of the dominant wavelength, optical bandwidth, and optical path difference of 
the interfering signals. Changes in the reflection spectrum caused by a periodic stretching and 
compression of the grating were experimentally measured using an unbalanced Michelson interferometer, 
a Michelson interferometer with a non-zero optical path difference. The interferometer’s sensitivity to 
changes in dominant wavelength of the interfering beams was measured as a function of interferometer 
unbalance and was compared to theoretical predictions. The theoretical analysis enables the user to 
determine the optimum performance for an unbalanced interferometer. 


Introduction 

Fiber optic Bragg gratings (FOBG) have become a useful tool for sensing applications. 1 ' 2 They have 
been successfully used in the validation of the structural integrity of bridges and aircraft components, 3-5 
temperature and pressure measurements in various applications, 6-8 and detection of leaks. 9 In addition to 
their use in static and quasi-static situations, dynamic changes in various system parameters have been 
measured using an FOBG. 10,11 In each of these applications, changes in a parameter of the system under 
observation such as strain or temperature causes a corresponding change in the wavelength of light 
reflected by the grating. A suitable detector detects the changes in the wavelength and generates 
corresponding signals that are further analyzed by signal processing electronics. 

This paper provides a theoretical foundation for analysis of optical signals reflected by FOBG’s under 
dynamic excitation. In particular, the performance of an unbalanced interferometer used as a spectrometer 
is analyzed, and the optical path difference producing maximal sensitivity to changes in the wavelength 
reflected by an FOBG is derived. An unequal optical pathlength interferometer was chosen because of its 
high sensitivity 12 as well as its compact size and low weight. These features make unbalanced 
interferometers especially attractive for in-flight health monitoring of aerospace vehicles and their 
components. 

The paper also provides a detail derivation of equations used to optimized the design of an 
unbalanced interferometer for the analysis of spectrally encoded signals. The equations are derived 
assuming the signals have a Gaussian power spectral density, and are consistent with a previously 
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reported analysis (Ref. 13). An experiment was devised to verify the theoretical dependence of the optical 
pathlength deference on the detected intensity. A preliminary version of this experiment has been reported 
elsewhere. 14 The experimental data are compared to the predicted results and possible causes of the small 
discrepancy between the two are discussed. 


Interference of Two Optical Beams in an Interferometer 

The light intensity produced by the interference of two electromagnetic beams depends on the 
intensities of the beams and their optical path difference, as well as on the dominant wavelength and 
bandwidth of each of the beams. Consider a beam entering an interferometer and having a dominant 
wavelength Ao and full-width-at-half-maximum optical bandwidth AA. After being equally divided by a 
beam splitter, the two component beams, each of intensity 7o, traverse two different paths, of length li and 
before being recombined at the detector. The optical path difference (unbalance) of the interferometer 
is nM = ri(li-l 2 ). The beams exiting the interferometer interfere at a detector with the intensity 
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where y(nAl) is the fringe visibility or fringe contrast reduction factor due to the spread in wavelength of 
each of the beams. We assume the spectral power density I(v) of each beam is Gaussian 14 and write the 
expression forlfv) using conventional notations 15 
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where Vo = c/2 o, c is the speed of light, and 
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The fringe visibility factor is the modulus of the complex degree of coherence, i.e., the Fourier 
transformation of the spectral power density function, 15 
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Using eq. (3) in eq. (4) we obtain 
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and substituting ct= nAl, the fringe visibility factor can be converted from the time domain to the spatial 
domain. 
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The intensity of the interference pattern at the detector then becomes 
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Figures 1 and 2 show I/I 0 as a function of the optical path difference for the ranges from 0 mm to 5 mm 
and 2.10 mm to 2.15 mm, respectively for Ao = 1300 nm and AX= 0.3 nm. These values are typical of the 
reflection spectrum of a generic optical communication wavelength FOBG. Note that Fig. 1 contains 
oscillations that are too fast to resolve. Those oscillations of a cosine nature are clearly seen in Fig. 2. 

As the parameters of the interfering beams change, the intensity of eq. (7) recorded by the detector 
changes as 
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If the optical path difference of the interferometer and the optical bandwidth of the beams remain constant 
while the dominant wavelength changes, the coefficient 
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gives the sensitivity of the intensity recorded by the detector to changes in dominant wavelength of the 
beams. The first term in square brackets in eq. (9) is the derivative of the rapid oscillations in Fig. 2 and 
the second term is the derivative of the much more slowly varying overall modulation in Fig. 1. Figures 3 
and 4 show the wavelength sensitivity function = 3(/// 0 )/9 Aq as a function of the optical path 


difference nAl for the values of Ao and AA used in Figs. 1 and 2. Since the second term in square brackets 
in eq. (9) is much smaller than the first term, the wavelength sensitivity function may be approximated by 
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If we had assumed instead that the optical path difference and dominant wavelength remained constant 
while the optical bandwidth varied, a similar expression for the bandwidth sensitivity function 

=d(l/l 0 )/dAX could be derived. The sensitivity function in eq. (10) possesses both 

rapid oscillatory dependence on nAl and slower linear and exponential dependence. In order to more 
easily compare with the experimental results of Sec. 3, we consider the envelope of the sensitivity 
function, which we call the envelope sensitivity function, 
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The function E ^ is graphed as a function of nAl in Fig. 5 for Ao = 1300 nm and AX = 0.3 nm. It is clear 

from Fig. 5 that in order to be sensitive to dynamic changes in the dominant wavelength of the signal 
entering the interferometer, the interferometer must have a nonzero unbalance or optical path difference. 

Setting the derivative of E^nAI) with respect to the optical path difference to zero, the path 
difference 


A; V21n2 

nAl = — 

max n aA (12) 

produces the maximum sensitivity of a given interferometer to changes in the dominant wavelength of the 
component beams. For example, for Ao and AA as in Figs. 1 to 4, we obtain flA/ max = 2.1 1 1 mm. The 
dependence of E^nAl) and tfA/ max on the optical bandwidth is shown in Fig. 6 for =1300 nm. The figure 

has plots of three curves that describe E^nAJ) and locations of its maximum that correspond to /7A/ max 
for three different values of AA. It could be also observed from eq. (12) that for a given wavelength Xo a 
decrease in the optical bandwidth AA leads to an increase in the value of /7A/ max . Because the bandwidth 
of optical signal is related to its temporal coherence, highly coherent optical signals would require the 
interferometer, in order for E^inAl) to reach its maximum, to have a very large unbalance / 7 A/ max . 
Practical considerations, in this case, should be used for the design of an appropriate interferometer. 

Since E% 0 (nAl) in Figs. 5 and 6 is slowly varying near / 7 A/ max , it is not necessary to operate the 
interferometer at exactly this unbalance in order to achieve good sensitivity. Rather, one can operate the 
interferometer within a relatively wide range of unbalance about /7A/ max and still obtain near-maximal 
sensitivity. This range of unbalance depends on the amount of deviation from the maximum sensitivity 
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that can be tolerated. We define the fractional deviation of E^nAl) from /: max as 
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To determine the range of imbalance for a given F we rewrite eq. (1 1), for simplicity, as: 
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where Xq is the value of x at which the curve E^inAl) reaches its maximum E ] 
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The resultant Taylor series in eq. (16) is approximated by retaining the first two nonzero terms, and after 
some algebraic simplifications we arrive at 
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Thus, for x in the vicinity of xo, the shape of the curve that describes E\Ax) in Figs. 5 and 6 resembles a 
parabola. Substituting eq. (17) into eq. (14) and rearranging the terms, we obtain 


F = 


£(*o )-£(*) 

E ( x o) 


= 2a 2 



V 


1 

a4l 


f 

J 



Replacing x and xo by the original argument nAl and ftA/ max , respectively, we yield 

\nAl — nA/ max (18) 

for the range of allowable imbalances about nAl max . If for example a sensitivity of the interferometer to 
within 1% of maximum is desired, an unbalance to within ±10% of the optimal value is acceptable. 
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Experiment 


The apparatus of Fig. 7 was used to measure the sensitivity envelope E^inM) of an unbalanced 

Michelson interferometer. A narrow portion of the -30 nm wide spectrum emitted by a superluminescent 
diode with the central wavelength of 1300 nm was reflected by a commercially available 
telecommunications-grade FOBG and entered the interferometer consisting of a cube beam-splitter and 
two mirrors, one of which was mounted on a translation stage. At the beginning of the experiment, the 
mirrors were positioned at slightly different distances from the respective nearest faces of the beam 
splitter. The FOBG was affixed to a lead zirconate titanate piezoelectric element (PZT) and a 1 KHz 
periodic voltage applied to the PZT generated vibrations that stretched and compressed the grating. The 
dominant wavelength reflected by the FOBG changed as a function of time in response to the periodic 
stretching and compression. The resulting time dependent light intensity exiting the interferometer was 
recorded by a PIN diode and was sent to an electronic spectrum analyzer where the amplitude of the PIN 
diode current at the PZT driving frequency was recorded. 

The amplitude of the detected signal at the PZT frequency was measured for a number of different 
values of interferometer unbalance. Changes in the optical path difference were produced by moving the 
mirror mounted on the translation stage. Each movement consisted of two steps. First, the mirror was 
moved approximately 250 fjun from its previous position. Next, the mirror was moved about this new 
position several times with an increment of about 0.1 jim in order determine the maximum amplitude of 
the PIN diode current at the PZT driving frequency. This maximum occurs at the quadrature point of the 
interferometer, and is proportional to Ex^inM). The measured wavelength sensitivity function is shown in 

Fig. 8. It is in good qualitative agreement with Fig. 5. Two similarly-shaped maximum sensitivity humps 
in Fig. 8 separated by a point of zero sensitivity are due to the fact that, in the process of changing the 
optical path difference (nAl), the movable mirror translated through a point corresponding to a zero 
optical path difference. This point is denoted in Fig. 8 as Z#. An obvious difference between Fig. 5 and 
the data of Fig. 8, however, is that although Fig. 5 predicts a linear variation of the sensitivity function 
near the zero path difference point, the variation in Fig. 8 appears to be roughly quadratic. In addition, 
the experimental data is more rapidly varying in the vicinity of the relative maximum than is the 
theoretical prediction. In particular, an optical path length 10% greater or less than the optimal values in 
Fig. 8 produces a sensitivity from 7% to 8% less than maximum, whereas the decrease in the theoretical 
sensitivity curve of eq. (1 1) is only 1%. 

To analyze the data of Fig. 8, we assumed the spectrum reflected by the Bragg grating was Gaussian 
and the optical bandwidth of the reflected spectrum was independent of the dominant wavelength. This 
last assumption has been found to be accurate for static measurements, 16 and we presume it is also the 
case for dynamic excitations of 1 KHz. The dominant wavelength of the reflected spectrum of the 
unstretched grating was determined using an optical spectrum analyzer to be Zo = 1310.2 nm ±0.1 nm, 
and agreed with the value Zo = 1310.17 nm ± 0.005 nm measured at the time of fabrication by the grating 
manufacturer. The optical bandwidth measured by the manufacturer at that time was AX =0.192 mil ± 
0.002 mil. The value of nA / max corresponding to this optical bandwidth via eq. (12) is 3.351 mm. 

In Fig. 8, the distances between the zero sensitivity point Z 0 and points of maximum sensitivity to the 
left and to the right of that point are about 1950 jim ± 25 jim and 2100 jam ± 25 \xm, respectively. 

Because of the double pass optical configuration of the Michelson interferometer, these distances are 
equal to half the optical path difference required to achieve maximum sensitivity of the interferometer to 
changes in dominant wavelength. Thus, the corresponding values of the optimal path difference /?A/ max 
are about 3.9 mm ± 0.05 mm and 4.2 mm ± 0.05 mm and the corresponding values of optical bandwidth 
via eq. (12) are AX = 0. 165 nm ± 0.002 nm and 0.153 nm ± 0.002 nm, respectively. The average of these 
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two values of the optical bandwidth is 17% below that quoted by the manufacturer. A least-squares fit of 
the results of Fig. 8 to the functional form of eq. (1 1) was not attempted. 


Discussion 

The good qualitative agreement between the theoretical sensitivity function and the experimental data 
underscores the correctness of our approach. The relatively small difference between our value of the 
optical bandwidth obtained experimentally, along with our theoretical prediction of the maximum 
sensitivity (eq. (12)), and the value provided by the manufacturer illustrates the basic validity of our 
theory. There are, however, a number of possible sources of error that could have contributed to the 17% 
difference. In particular, the interferometer was not locked at the quadrature position. As a result, slow 
thermal drifts of the optical path difference or building vibrations could have introduced errors into the 
measurements. A solution to this would be construction of an unbalanced interferometer on a chip with 
integrated compensation. Also, in addition to assuming the reflection spectrum of the FOBG was 
Gaussian and the optical bandwidth was independent of dominant wavelength for dynamic excitations, we 
also ignored the variation in the -30 nm wide superluminescent diode spectrum over the ~0.2 nm wide 
bandwidth of the FOBG reflection spectrum. This effect, however, is expected to be small. 

The slowness of variation of the theoretical sensitivity curve in the vicinity of the relative maximum 
works to one’s advantage when designing an FOBG sensor system. The value of optimum optical 
difference nAl max obtained experimentally (Fig. 8) varies roughly from 16% to 25% of the unbalance 
77A/ max = 3.351 mm, the value predicted by eq. (12) using the FOBG manufacturer’s optical bandwidth 
measurement. As a result, one can use eq. (12) as a reasonable estimate for setting the interferometer 
unbalance since small deviations of the actual optimal unbalance point from the prediction of eq. (12) 
produce correspondingly small decreases in the sensitivity of the interferometer to changes in the 
dominant wavelength of the input signal. 
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Optica! path difference, nil, mm 


Figure 1 . — Intensity of two interfering optica! beams of 
dominant wavelength Xq ~ 1300 nm and FWHM optica! 
bandwidth AX = 0.3 nm as a function of optica! path 
difference. 
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Micrometer readout, microns 


Figure 8. — Measured power at the PZT driving fre- 
quency as a function of the micrometer readout 
of the movable arm of the interferometer. 
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